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Abstract. In this paper we summarize our recent results for low energy associated 
charmonium production cross sections, using 1) crossing symmetry, and 2) an explicit hadronic 
model. These predictions are of relevance to the planned charmonium and charmonium hybrid 
production experiment PANDA at GSI. 



1. Introduction 

The clear spectrum of states observed in the charmonium system and the straightforward 
identification of low-lying states with the predictions of simple cc potential models make this an 
ideal sector for the search for unusual "extra" states such as charmonium hybrids. The PANDA 
experiment pQ at GSI plans to exploit this possibility in a search for charmonia, charmonium 
hybrids, and other similar states in a pp annihilation experiment. 

Of course planning for this experiment would be greatly facilitated by experimental data on 
and theoretical estimates of these cross sections. Unfortunately, remarkably little is known at 
present. There are a few data points for the single reaction pp — » 7r° J/ip from the E760 |4j and 
E835 [6] experiments at Fermilab, which suggest a scale of roughly 100 pb near 3.5 GeV, but 
nothing is known about the cross sections to other charmonium states. Here we give results for 
two estimates of these cross sections, based on two different theoretical approaches. Where they 
allow common predictions, we find that these methods are reassuringly in agreement to within 
about a factor of two, and the explicit hadronic model shows a similar level of agreement with 
the data. The following discussion summarizes these predictions in more detail. 



2. Crossing 

This approach was developed by Lundborg et al. [2j. The close kinematic proximity of certain 
charmonium processes with identical underlying amplitudes, such as the production process 
pp — > 7r°J/ip and the experimentally well studied three body decay J/ip —> 7r°pp, suggests that 
knowledge of the decay can be used to estimate the numerical scale of the production cross 
section. This is illustrated by FigJTJ in which it is evident that an extrapolation from the decay 
Dalitz plot to the production regime may well be feasible. 

The differential decay rate and differential production cross section are the following functions 
(with an implicit polarization sum) of the same invariant amplitude M\ 
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Figure 1. Kinematically allowed regimes for the charmonium production reaction pp — > 7r° <//?/> 
and the related decay J/ip — ► TT°pp. 
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Here m is a generic meson with spin S m and ^ is a generic charmonium (or related) state. 
Ideally one can develop a detailed model of M from a study of the decay Dalitz plot, which 
can then be extrapolated to predict the charmonium production cross section. However a much 
simpler qualitative cross section estimate can be written if the Dalitz plot does not show much 
structure, which is the constant amplitude approximation. If we assume that is constant, 

it can be eliminated between the decay and cross section formulas, which gives a very simple 
cross section estimate in terms of the three-body decay rate: 
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where is the area of the decay Dalitz plot. 

This simple constant-amplitude prediction for the cross section for pp — ► 7r° 'J ftp is shown in 
Figj2j This is the single experimentally measured case; as we noted previously, there are a few 
data points for this cross section from the E760 [lj and E835 [5j[6] experiments at Fermilab, which 
are also shown in the figure. Evidently there is indeed rough agreement between theory and 
experiment, which suggests that extrapolation from data on generic decays of the form \E' — > mpp 
to the corresponding associated production cross section pp — * m^f is a useful approach. 

The corresponding results for a range of J ftp and ip' decays are given in Ref.pj. One 
concern in applying this approach is that some decay processes may be dominated by strong 
N* resonance bands, in which case extrapolation to the production cross section would be 
inaccurate. An example is J/ip — * rjpp, which presumably has large contributions from the 
iV*(1535). Extrapolation of this and other decay processes to production cross section should 
generally be accompanied by careful studies of the decay Dalitz plot which allow N* resonance 
contributions. 
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Figure 2. Constant-amplitude prediction for o~(pp — > TT°J/ip). The data points are Fermilab 
measurements from E760 [4J (filled) and E835 [5] (open). There are additional experimental 
results for this reaction from E835 [6] which have not yet been published as a physical cross 
section. 



3. Hadronic Model 

An attractively simple dynamical model which we have recently developed [7j assumes that 
the associated production of a charmonium or charmonium hybrid state ^ with a 7r° in pp 
annihilation takes place through pion emission from an incoming p or p line. The process 
pp — ► vr ^ at leading order then involves the two Feynman diagrams of Figj3j in which the only 
a priori unknown quantity is the ^ — pp vertex. 




Figure 3. Feynman diagrams assumed in this model of the generic reaction pp — ► ir 

The use of a similar model, with a PCAC axial vector coupling for the ppir vertex, was 
earlier suggested by Gaillard, Maiani and Petronzio [3J. 

We have used this model (with a simpler pseudoscalar ppn vertex) to calculate the cross 
sections for pp — > ir°fy for a much wider range of light charmonia, ^ = r) c , J/ip, XOi Xi an d ip'- 
Numerical results for these cross sections require the strength of the coupling ^ — pp, which 
we take from the measured annihilation widths to pp. (We do not consider the X2 because of 
the more complicated vertex, and other charmonium states do not have known pp couplings.) 
Analytical results for all these total and differential cross sections and the numerical results 
shown here are given in Ref.[7]. 

Our results for these cross sections are shown in FigHJ First, note that the single measured 
case of pp — > 7r J/ip is indeed in rough (factor of two) agreement with the predictions of this 
dynamical model. This suggests that the simple Feynman diagrams of Fig|3] may well give a 
reasonably accurate description of these processes, which is very encouraging. Application of 
the model to other charmonium and charmonium hybrid states only requires knowledge of the 
coupling ^ — pp. 
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Figure 4. Predicted unpolarized total cross sections for the processes pp — > vr ^, where 
$ = r/ c , J/ip, Xo,Xi an d The data points are the E760 (filled) and E835 (open) experimental 
results for ^ = J/ip. 



Figure 5. The corresponding angular distribution for the reaction pp — > tt rj c , for E cm = 
3.2 - 5.0 GeV in steps of 0.2 GeV. 

The most remarkable predictions of the model are the very large cross sections for the 
associated production of the rj c and xo states relative to the J ftp (see FigJl]). In particular, the 
cross section for pp — > ir r] c is predicted to be about a factor of 30 larger than the pp — > -K G J/ijj 
cross section in the range of yfs relevant to PANDA. The reason for the large r] c and xo couplings 
to pp is not known, but it may be that pp annihilation through gg final states is significantly 
larger than the ggg states that couple to the J /ip. This suggests that production of some hybrid 
states may be enhanced at PANDA over cc, since gg can couple to the hybrid basis state ccg 
more easily than to cc itself |8j. 

This model also gives interesting predictions for the angular distributions of charmonia 
produced in pp associated with a 7r°. As an example, the cm. frame angular distribution 
for the largest cross section, pp — > 7r°r/ c , is shown in Figj5j this angular distribution has a very 
characteristic node at 9 = tt/2, due to a t «-> u antisymmetry in this model. Since the angular 
distributions we find for pp — > vr ^ are sensitive to the quantum numbers of the charmonium 



state they could be used as tests of the accuracy of this simple hadronic model. 
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